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ABSTRACT 

The symmetries of the DNA double helix require a new term in its linear response 
to stress: the coupling between twist and stretch. Recent experiments with torsionally- 
constrained single molecules give the first direct measurement of this important material 
parameter. We extract its value from a recent experiment of Strick, et al. and find rough 
agreement with an independent experimental estimate recently given by Marko. We also 
present a very simple microscopic theory predicting a value comparable to the one observed. 

INTRODUCTION 

In this paper we will study the response of DNA to mechanical stress using the methods 
of classical elasticity theory|jl|. While many elements of DNA function require detailed 
understanding of specific chemical bonds (for example the binding of small ligands), still 
others are quite nonspecific. Moreover, since the helix repeat distance of /q ~ 3.4nm involves 
dozens of atoms, it is reasonable to hope that this length-scale regime would be long enough 
so that the cooperative response of many atoms would justify the use of a continuum, classical 
theory, yet short enough that the spatial structure of DNA matters. 

Since various important biological processes involve length scales comparable to Iq (no- 
tably the winding of DNA onto histones), the details of this elasticity theory should prove 
important. Yet, until recently, little was known about the relevant elastic constants. Exten- 
sive experimental work yielded fair agreement on the values of the bend and twist persistence 
lengths, though the former was plagued with uncertainties due to the polyelectrolyte charac- 
ter of DNA . A simple model of DNA as a circular elastic rod gives a reasonable account 
of many features of its long-scale behavior, for example supercoiling 0. 

Recently, techniques of micromanipulation via optical tweezers and magnetic beads have 
yielded improved values for the bend stiffness from the phenomenon of thermally-induced 
entropic elasticity @|^[^, as well as a direct measurement of a third elastic constant, the 
stretch modulus00. Significantly, the relation between bending stiffness, stretch modulus, 
and the diameter of DNA turned out to be roughly as predicted from the classical theory of 
beam elasticity |3I[@]0) supporting the expectations mentioned above. 

Still missing, however, has been any direct measurement of the elastic constants reflecting 
the chiral {i.e. helical) character of DNA. One such constant, a twist-bend coupling, was 



investigated by Marko and Siggia []10| , but no direct experimental measurement has yet been 
devised. We will introduce a new chiral coupling, the twist-stretch energy. Electrostatic 
effects do not complicate the analysis of this coupling. We will explain why our term is 
needed, extract its value from the experiment of Strick, et al. [|TI|], and compare it to the 
prediction of a microscopic model to see that its magnitude is in line with the expectations 
of classical elasticity theory. J. Marko has independently introduced the same coupling and 



estimated its value from different experiments] 12]; our values are in rough agreement. 



EXPERIMENT 



DNA differs from simpler polymers in that it can resist twisting, but it is not easy to 
measure this effect directly due to the difficulty of applying external torques to a single 
molecule. The first single-molecule stretching experiments constrained only the locations of 
the two ends of the DNA strand. The unique feature of the experiment of Strick et al. was 
its added ability to constrain the orientation of each end of the molecule. 

We will study Fig. 3a of ref. In this experiment, a constant force of 8pN was applied 



to the molecule, and the end-to-end length Ztot was monitored as the terminal end was rotated 
through ALk turns from its relaxed state (which has Lk^ turns). In this way the helix could 
be over- or undertwisted by as much as ±10%. Over this range of imposed linkage z^^t was 
found to be a linear function of a: 

e = e^=o- 0.15(7, (1) 

where a = ALk/Lko and e = (ztot/ztot,o) ~1- Thus a is the fractional excess link, and e is the 
extension relative to the relaxed state. Eq. ([^) is the experimentally observed twist-stretch 
coupling. 

THEORY 

Phenomenological Model 

A straight rod under tension and torque will stretch and twist. We can describe it using 
the following reduced elastic free energy per equihbrium length Ztot,o of the rod: 
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/i((T, e) = , ' ^ = ^ Ca' + Be' + 2Dea . (2) 



The twist persistence length is (7 ~ 75 nm 0, while the helix parameter ujq = 27r//o = 
1.85/nm. We will take B ^ llOOpN/uj^kBT ^ 78nm In the experiment under study, 
there is an applied reduced force r = SpN/fc^T ^ 1.95/nm. For a circular beam made 
of isotropic material, the cross-term D is absent ^ because twisting is odd under spatial 
inversion while stretching is even. For a helical beam, however, we must expect to find this 
term. 

Setting r = dfi/de\a, we find 

e = e^=o - iD/B)a . (3) 

Comparing to Eq. (|I|), we obtain the desired result: D = 12 nm. To compare this to Marko's 
analysis, we note that his dimensionless g equals our Dujq, so that we get g = 22. The rough 
agreement with Marko's result g = 35 |0 indicates that the data show a real material 
parameter of DNA and not some artifact. We do not expect exact agreement, since Marko 
studied the nonlinear overstretching transition of [0]l3; our value came from the linear 
regime of small strains. 

Microscopic Model 



To gain further confidence in our result, we will now see how the expected twist-stretch 
coupling emerges from a simple elastic model for DNA. Fig. |1| reviews the relevant geometric 



(a) (b) 




E2 

(c) 

Figure 1: (a) Schematic representation of a cross section of DNA showing its intersections (Bi 
and B2) with the phosphate backbones, its intersection (hne A1PA2) with the parameterizing 
hehcal ribbon, and its hehcal center C. (b) A section of ribbon showing oppositely oriented 
edges. At each point along the ribbon's center (dotted curve) , there is a triad of orthonormal 
vectors, Ei, E2, and E3. E3 is parallel to the center line, E2 points to one of the edges, and 
El is perpendicular to the ribbon, (c) Representation of the helical ribbon. 



properties of DNA. Base pairs are connected by inequivalent sugar-phosphate backbones 
whose twisting pattern defines the major and minor grooves. These backbones intersect a 
given cross section of DNA at two points. Lines parallel to the line connecting these two 
points sweep out ribbon-like surfaces as the DNA twists. We will parameterize the DNA 
structure by one of these ribbons, whose center is a distance Tq from the central helical axis. 
Any choice of Tq is acceptable. We will see, however, that a small value of Tq can explain the 
measured value of D. 

Our ribbon is described by the triad of unit vectors Ei, E3, and E2 = E3 x Ei, where E3 
is the unit tangent vector to the center of the ribbon and E2 points from the center of the 
ribbon to one of its edges. The triad varies as we move along the arc length s of the ribbon. 
The motion is described by 

dEj 

— —^ijk^j^k- (4) 

The parameter s labels each point along the central axis of the ribbon in its unstressed state 
and runs from to L. The actual arc length along the distorted central axis of the ribbon 
will not be ds but rather [1 -|- e2{s)]ds where €2 is the intrinsic strain. Therefore, the total 
length for constant €2 is L = (1 -|- €2)!/. The intrinsic strain allows the spacing between 
successive phosphate groups to change. 

The edges of our ribbon, like the two sugar-phosphate backbones in DNA, are distin- 
guishable and point in opposite directions. This symmetry can be incorporated into our 
microscopic model by considering a rotation of 180° about the vector Ei followed by s — ^ — s. 
Under this transformation E2 and E3 change sign, but s derivatives of these vectors do not. 
Also, El does not change sign, but dEi/ds does. Therefore, the free energy should remain 



unchanged upon changing the sign of f2i but not of and i^sIT^. The most general reduced 
free energy per length of ribbon relative to that of the flat unstretched ribbon up to second 
order in Qi and 62 is 



Idna = ^[A'Q,' + A{Q2-^2of + CiQs-^3oY + BuJo'el + 

62 does not change sign under s — > — s, and so it can appear in combination with ^2 and Vt^ 
in Eq. (1). This model is the simplest semi-microscopic model that incorporates all of the 
symmetries of DNA. It is an expansion to harmonic order in first-order derivatives of the 
vectors Ej [e.g. of ^3 = — Ei ■ Thus, it is a model with quantitative predictive power 

so long as the VLi are slow on a scale set by the distance a = 0.6 nm between successive 
phosphate groups, {i.e. so long as fijO <^ 1). In the ground state Q3 = Q^o ~ cuq = 1.85/nm, 
so that a ^ 1.1 is not small. This implies that higher derivative terms {e.g. {^^y, 
etc.) are needed for a quantitative theory. Nevertheless, our simple semi-microscopic model 
captures the essential symmetry of the DNA structure and allows us to address questions 
like the nature of the twist-stretch coupling. 

One can easily show that the center of the ribbon describes a helix in the ground state of 
Eq. d^). We will assume that D, G, and K can be made small by an appropriate choice of 
tq. Then, to keep the model as simple as possible, we will simply set these parameters equal 
to zero for this choice of tq. We parameterize the helical ribbon using three angles ip, 7, and 
0: 

E3 = sin 7Z + cos 7(/) 

E2 = cos '?/'(— sin 70 + cos 72) — sin '?/'/) 

El = — sin?/'(— sin70 + cos7z) — cos'j/'p, (6) 

where p and are cylindrical unit vectors spinning at frequency 0. If is a constant and 
■0 = 0, then the ribbon wraps around a cylinder of length Ztot and radius r. 

cos 7 

Ztot = L{^ + £2) sin 7 and r = . (7) 

00 

In its ground state the helix has the following properties: Qi = Qiq = 0, ^2 = ^20 = 1^0 cos 70, 
^3 = ^30 = Co'osin7o, = cuq = \J^2o^ + ^30^5 and -0 = -00 = 0. Also, the ground state 
length and radius of the molecule are, respectively, 2tot,o = -Z^ sin 70 and ro = cos7o/ci^o- 

We now consider deviations in the ground state length Ztot and twist rate of the 
molecule. Since the total twist is 0L, the excess twist is aL where a = {(j) — ojq) / ojq = Scp/uQ. 
Using Eq. (|^) we find that changes in length are produced both by intrinsic strain 62 and by 
changes in 7: 

e = — 1 = ei + 62 where ei = cot 70^7. (8) 

^tot,0 

The energy of harmonic deviations from equilibrium are obtained by expanding J'dna to 
second order in Sj{s), Sip{s), and 50(s). The ground state is a periodic helix implying these 
variables can be expressed in terms of Fourier modes in different Brillouin zones defined by 
ujq. Rotations of the helix about space-fixed axes x, y, and z are described by the variables 
66x, S6y, and 66z. The variables ^7, 6tp, and 50 can be expressed in terms of these variables: 



^7 = cos cuqs 59 X + sin uoqs 59 y 

5ip = (sin uqs 59x — cos uqs 59y) 

COS70 

(50 = tan7o(sinc<Jos (56*^ — coscuqs (56*^) + 56*^. (9) 

Thus, variations of 59^ and 59y in the first BZ give rise to variations in ^7 and 5il) in the 
second BZ. A complete long-wavelength theory can, therefore, be expressed in terms of the 
first BZ components of 59x-i 59y, ^7, 5ip, and 5(j) (whose 1st BZ component is equal to that 
of 59,). 

Using the relations for the Qi obtained from Eq. (^) and Eq. (|^), we find 

fl2 — ^20 = cos 7o50 — ujo sin 70^7 

n^-ilso = sin7o50 + u;ocos7o57. (10) 

Then, using these expressions in Jdna, integrating out tp, and remembering that length 
along the pitch axis is a factor of sin 70 less than the ribbon length, we find that the effective 
reduced free energy per unit length of pitch axis is 



where 



f = fB + fTS, (11) 



fs = T^iA' + Asin^o + Ccos^o)(^^ + = \a{91 + ^J) (12) 
4sm7o / 



is the bending energy [P!B| and /t5 is the twist-stretch energy defined by 



hs = TT^^iCaaCr^ + 5,,eiei' + Be2^ + 2D,,^eia], (13) 
/ sm 7o 

where a,. = (Acos2 7o + Csin2 7o), B,^,^ = {A'^^^ + C sin^ -fo) , and D,^^ = {C-A)sm^-ro. 
The twist-stretch energy can be expressed in terms of the total strain by setting ei = e — €2 
and integrating over €2- The result is that Jts has the same form as Eq. (^ with macroscopic 
elastic constants B, C, and D expressed in terms of our microscopic parameters A, B, and 
C. The stretch moduli associated with ei and €2 add in parallel to yield a total stretch 
modulus B = h 4)~^- The twist-stretch coupling is 

D= ,^ ^ (C-A)sin7o- (14) 

These expressions for the twist modulus stretch modulus 5, and twist-stretch modulus 
D are valid for arbitrary values of the helix offset tq. 

Thus, a description of DNA in terms of a helical ribbon with an axis offset from the 
central helical axis generates a twist-stretch coupling even if the bare twist-stretch coupling 
{D in Eq. (|)) is zero. We can estimate the offset tq necessary to produce the measured D 
assuming it arises entirely from Eq. (p!^ . If we assume A' ^ A and ujqTq -C 1, then B ^ B = 
75nm, C = 78nm, and A ^A = 40nm, and we find D ^ {(jJoro)'^(B/A)(C - A) and 
(ujQroY ^ 0.176 or tq ~ 0.23nm[|l5l. Corrections to this estimate are of order (tuoro)'^ ~ 0.03. 



This result for the twist-stretch couphng D in the hmit of a small helix offset ro was found 
previously in 

CONCLUSION 



We have modeled DNA as a thin helical ribbon and presented a complete long-wavelength 
theory which includes energy costs due to bending, stretching, and twisting the DNA molecule. 
Using this theory, we were able to relate semi-microscopic elastic constants to the experi- 
mentally measured macroscopic elastic constants. We have also calculated the coefficient of 
the twist-stretch coupling and compared it to torsionally constrained DNA stretching ex- 
periments. We have found that the experimental value of the twist-stretch coupling gives a 
value ro for the radius of the helical ribbon in rough agreement with the elastic center offset 



from the helix axis obtained from crystallographic data[|T5 
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